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Abstract. A planar slab of negative index material works as a superlens with sub- 
diffraction- limited imaging resolution, since propagating waves are focused and, moreover, 
evanescent waves are reconstructed in the image plane. Here, we demonstrate a superlens 
for electric evanescent fields with low losses using perovskites in the mid-infrared regime. 
The combination of near-field microscopy with a tunable free-electron laser allows us to 
address precisely the polariton modes, which are critical for super-resolution imaging. 
We spectrally study the lateral and vertical distributions of evanescent waves around the 

^-H image plane of such a lens, and achieve imaging resolution of A/ 14 at the superlensing 

JL/ wavelength. Interestingly, at certain distances between the probe and sample surface, we 

, observe a maximum of these evanescent fields. Comparisons with numerical simulations 

indicate that this maximum originates from an enhanced coupling between probe and 
^ object, which might be applicable for multifunctional circuits, infrared spectroscopy, and 

thermal sensors. 
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Introduction 

In 1968 Veselago theoretically discussed a material with a negative refractive index, as 
well as its extraordinary responses to electromagnetic waves including negative refraction, 
reversed Doppler and Cerenkov effects and focusing with a planar lens [Ij. Pendry revis- 
ited this idea 30 years later, showing that with such a planar lens one can create an image 
with a resolution beyond the conventional diffraction limit [2]. This opened the rapidly 
growing field of metamaterials, which present opportunities for new and remarkable appli- 
cations in electromagnetics Oil] ranging from electrostatics [5^ 6J via radio frequency 
microwave [8], and Terahertz [9] to the infrared [lOj and optical regime [HI [121 [El E]. 
The experimental proofs of negative refraction and sub-diffraction-limited resolution by a 
negative-index planar lens were first reported in the GHz range [TSl [16]. Subsequently, 
it was shown for visible [17j and infrared (IR) [18j wavelengths that the evanescent field 
information of transverse-magnetic (TM) waves can be recovered using a medium with 
negative permittivity (e) only. Such a lens is normally called a superlens, and substantial 
progress has been achieved in this rapidly developing area [19j . 

A superlens is preferably realized by low-loss materials [2l [T9] , among which perovskite 
oxides are good candidates. In the mid-infrared (mid-IR) range, perovskite oxides show 
phonon resonances, which are located at slightly different wavelengths for related materials 
such as e.g. bismuth ferrite (BiFeOa) and strontium titanate (SrTiOa) [20, |2T]. On the 
high-frequency side of these phonon resonances, perovskites exhibit negative permittivities, 
which are suitable for superlens structures [21 [19]. At these wavelengths, the intrinsic 
absorption of light is small [20l [21] . In addition, epitaxially grown oxides exhibit highly 
crystalline interfaces resulting in low scattering. Both the small intrinsic and low scattering 
losses of perovskites could significantly improve the imaging resolution of superlenses. 

Perovskite oxides exhibit many intriguing properties such as colossal magnet or esistance 
[22], ferroelectricity [23j, superconductivity [24J, and spin-dependent transport [25j, which 
lead to numerous applications such as non- volatile memories microsensors and mi- 
croactuators [27], as well as in nanoelectronics [28j. Piezoelectricity and ferroelectricity 
allow for the manipulation of geometrical dimension, electric polarization, and dielectric 
properties by external electric fields [29j. Multiferroic perovskites, such as BiFeOa, exhibit 
ferroelectricity as well as magnetoelectric coupling [29], thus providing a pathway for addi- 
tional degrees of tunability. Moreover, perovskites with matching lattice constants can be 
grown epitaxially on top of each other |30^ 131] . which allows one to combine their proper- 
ties in multifunctional heterostructures [^ ^] . These unique properties of perovskites 
may inspire new functionalities of metamaterial structures such as tunability by external 
fields and novel designs for multifunctional circuits. 

In this paper, we study a new type of superlens for electric fields in the mid-IR based 
on the perovskite oxides BiFeOa and SrTiOa.We investigate the evanescent waves in the 
image plane of perovskite superlenses by means of scattering-type near-field infrared mi- 
croscopy (s-NSIM) [311 [35]. The combination with a free-electron laser (PEL), which is 
precisely tunable in the wavelength (A) regime from 4 to 250 /xm, enables us to address 
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the polariton modes and study their lateral, vertical and spectral distribution. Such de- 
tailed characterizations are important to understand superlensing effect, but they were not 
comprehensively conducted in previous work [TT^ "18]. At certain wavelengths we observe 
enhanced evanescent fields in the image plane as well as a resolution beyond the classical 
diffraction limit. Moreover, we find that the evanescent fields show maxima at a certain 
distance between the probe and the sample surface, for which the vertical position depends 
on the wavelength. Comparisons with two-dimensional numerical simulations indicate that 
a superlens-enhanced coupling between probe and object causes this effect, which might 
be utilizable for controlling polariton propagation. 

Results 

Perovskite-based superlens. The original superlens proposed by Pendry is a single slab 
of a material with ^ = — 1 surrounded by air {eair = +1) [2J. At the two interfaces of the 
layer coupled surface polariton modes are excited, which amplify the evanescent electric 
fields arising from an object and transform them to the opposite side of the lens [2J. The 
oxide-based superlenses in our study consist of layers of matched perovskites (see Fig. la). 
We structure a 50 nm thick film of metallic perovskite strontium ruthenate (SrRuOa) on a 
SrTiOs substrate by photo-lithography, acting as superlens objects with a lateral size down 
to 3 fim. In order to ensure a fixed distance between objects and the slab we add a layer 
A (thickness d, £a — ^air — +1) on top of the objects before we grow the actual superlens 
layer with = — 1 (thickness 2d). Since layer A has a slightly different permittivity 
than air this two-layer superlens is asymmetric. Hence, we additionally study a superlens 
with symmetric design after adding an additional layer of material A (thickness d) on top. 
Superlensing is expected for such systems, when the real parts of the permittivities ^e{£A) 
and ^e{£B) have the same absolute values and opposite signs 3fte(5^) = —^e{£A) < 0. 
Assuming the dielectric constants of the layers are similar to the ones described in literature 
[201 El] (see Supplementary Fig. SI), this superlensing condition is fulfilled at A = 13.9 /xm 
with the layers A and B being BiFeOa and SrTiOa, respectively (Fig. lb,c). Many other 
similarly suitable pairs of perovskites, such as PbZra;Tii_a;03 and SrTiOa, BiFeOa and 
PbZr^Tii_^03, BaTiOa and SrTiOa, or BiFeOa and BaTiOa, can be found at different 
wavelengths, implying the flexibility of perovskite superlenses in terms of the operation 
wavelength. Please note that other groups of dielectrics such as e.g. fluorides |36l ET] and 
simple oxides [38] can form similar pairs for superlenses, which altogether cover an even 
larger wavelength range of operation. 

We examine the enhanced evanescent fields on the image side of the lens with s-NSIM 
in combination with a free electron laser (FEL) (see Fig. la and Methods), which allows 
for polariton-enhanced s-NSIM Surface polariton modes can be excited at the 

interface of two materials which have permittivities of opposite signs '^2]. In the 
studied wavelength range from 13.5 fim to 18.5 /xm, the layers A and B of the perovskite- 
based superlenses show the following permittivities (see Figs. lb,c) and corresponding 
polariton modes: When A > 16.5 fim we find '^{eb) < ^(^a) < and polariton modes are 
excited at the entire sample surface, being the interface between the toplayer and air. For 
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A < 16.5 /xm it is ^{sb) < < ^(sa) and polariton modes are created at both interfaces 
of layer B. As layer B is only 400 nm thick, these two modes are coupled, but are not 
necessarily localized. We observe an enhanced near-field signal, whenever a polariton mode 
is present at the position of the probe. However, only around the superlensing wavelength 
{eA — —^Bi) the polartion modes are strongly localized and, hence, create a sub-diffraction 
hmited image of the objects on the opposite side of the lens. Please note that with 'near- 
field image', we refer to the image obtained by near- field microscopy. The contrast in 
these images arises from local sample properties and thus different coupling to the probe, 
allowing us to distinguish the different objects on the opposite side of the superlens with 
a sub-diffraction-hmited resolution (see Methods) [43j. 

Mid-IR near-field imaging. When placing the probe at a distance of about 30 nm to 
the sample, which is scanned relatively to the probe position, we obtain the near-field 
response as well as the topography of the sample. The corresponding results are depicted 
in Fig. 2 for similar object distributions with lateral sizes from 3x3 fiw? to 8x8 jiw? and 
three different samples, namely (a) structured SrRuOa objects on a SrTiOa substrate, (b) 
a symmetric superlens and (c) an asymmetric superlens without the top BiFeOa layer. 

The SrRuOs objects without a lens (Fig. 2a) show a near-field response, which depends 
on the polarization and the incident angle of the impinging light with respect to the objects. 
The fc-vector direction of the incident TM-polarized light is shown in the sketch and in 
the topography image in (a) and is the same for all scans depicted in Fig. 2. Over a wide 
wavelength range in the IR, the permittivity of SrRuOa is much smaller than zero, resulting 
in an excellent, metal-like scattering behavior. Consequently, the distance dependence of 
the near-field signal on SrRuOa exhibits the same exponentially decaying character of the 
evanescent waves for all wavelengths (see middle panel of Fig. 2a). The metallic SrRuOa 
structures show enhanced field intensities when their geometries match the wavelength, an 
effect which is well known in plasmonics [44J. Consequently, certain parts of the objects 
appear bright in the second- and third-harmonic near-field signals NF2Q and NF^q (see 
Methods). In general NF2Q is about 3 times larger than NF^q, with an interference-like 
background signal. In order to ensure pure near-field detection in the following, we will 
depict third-harmonic signals only. 

For the symmetric superlens (Fig. 2b), at A = 17.5 /xm we observe a near- field signal 
due to a non-localized polariton mode at the sample surface, showing no clear contrast 
between areas with and without SrRuOa objects. Close to the superlensing wavelength, at 
around 14.6 /im, the evanescent fields are highly localized, resulting in a strong contrast 
reconstructing the SrRuOa structures. Please note that even though the NSIM probe scans 
more than 800 nm higher above the SrRuOa objects compared with the bare SrRuOa case, 
the signal is about 2 times stronger because the coupled polariton modes at the interfaces 
enhance the near-field signals arising from the objects. 

On the asymmetric superlens (Fig. 2c) we observe a similar, slightly blue-shifted response 
compared to the symmetric superlens, with non-localized polariton modes around 17.3 /im 
and superlensing with maximum contrast at around 14.1 fim. The signal on the asymmetric 
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super lens is 2 times stronger than for the symmetric super lens, as the near-field probe is 
placed by 200 nm closer to the top SrTiOs surface at which the polariton mode is excited. 

With both superlenses we clearly resolve the buried SrRuOs objects at the corresponding 
superlensing wavelengths. In the following, we discuss in detail the spectral response of 
both samples over a broad wavelength range. 

Spectral response of the near-field signals. Our s-SNIM setup allows for studying 
polariton modes with different characteristics in the wavelength range from 13.5 jj^m to 
18.4 fim. For the symmetric superlens (Fig. 3), we observe an enhanced localized signal for 
A from 14 to 15.5 /im, close to the superlensing condition. The highest contrast is observed 
at 14.6 /im and the signal drops to zero for wavelength 1 /im smaller or larger than that. 
Note that areas without objects show a small negative signal in the superlensing regime, 
possibly due to evanescent waves on the substrate superlens interface with opposite phase. 
At longer wavelengths both layers, BiFeOs and SrTiOs, show negative permittivities and 
a single polariton mode is excited at the topmost sample surface leading to an enhanced 
near-field signal on the entire sample area for A > 17 /im. Even though this signal is up 
to 8 times higher than the superlensed signal, it does not show any image contrast except 
for some artifacts at the topography edges caused by scattering and errors in the distance 
control. 

For the asymmetric superlenses (Fig. 4), we observe a localized near-field signal around 
the superlensing wavelength for 13.8 /im < A < 15 /im with a clear object-related contrast 
at 14.1 /im. In addition, the top-most layer SrTiOs of this lens supports a non-localized 
polariton mode for 14 /im < A < 15.5 /im. In contrast to the symmetric case, in which this 
mode is spectrally well separated from the superlensing wavelength, for the asymmetric 
lens the wavelength regimes of superlensing and non- localized polariton overlap. Interest- 
ingly both modes respond with opposite phases resulting in a destructive superposition on 
the objects, which appear dark in the range from 14.5 /im to 15.2 /im. With increasing 
wavelength the non-localized mode dominates and the images are blurred. 

The smallest structures resolved on the superlenses are about 3x3 /xm^ in size and about 
1 /im apart corresponding to a relative resolution of about A/ 14. This sub-diffraction lim- 
ited image is created by the superlensing effect. Compared to the image of bare SrRuOs 
objects the superlensed image appears homogeneously bright as the incident light is scat- 
tered by the probe and from there directed towards the object (see Methods). In general, 
the imaging with s-NSIM is influenced by the probe and its vertical position as it is in 
particular observed in the cross sections of Fig. 5, which will be discussed in the following 
paragraphs. 

Mid-IR near-field cross sections and spectroscopy. We focus on the asymmetric 
superlens to study the vertical distributions of evanescent waves. Similar measurements 
for the symmetric superlens can be found in the Supplementary Fig. S2. Fig. 5a depicts 
normalized vertical cross sections on the image side of the sample for different wavelengths 
showing the near-field signal as a function of the probe-sample distance z and the position 
of the sample. The topography of the sample is reflected by the dark area on the bottom 
of the pictures with a 4 /im-wide SrRuOa object in the center. For A > 14.8 /im the probe 
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excites propagating polariton modes close to the SrTiOa surface. Around the superlensing 
wavelength, for A = 14.5 to 13.5 /im, the evanescent field is locahzed on the object with 
decreasing signal and contrast for shorter wavelengths. For the latter wavelength regime, we 
observe an intriguing phenomenon: a maximum in the evanescent field appears at a certain 
distance zq between tip and sample surface. Moreover, zq increases with smaller wavelength 
and appears at distances of up to 150 nm from the sample surface at A = 13.5 fim. At 
first glance, this effect seems to be unexpected because the phonon polariton mode is a 
confined surface mode, which exponentially decays from the interface. To make sure that 
this effect is not an artifact in our experimental setup, we compare these results with 
numerical simulations as discussed in the following. 

Fig. 5b shows the simulated data for a superlens consisting of the same constituents and 
geometry as in our experiments (see Methods) . In these simulations we clearly observe the 
same effect as in Fig. 5a, that is, the field maximum gradually shifts away from the sample 
surface when the wavelength decreases. This maximum has an asymmetric shape, locating 
on the right-hand side of the structure, which is illuminated from the left with an incident 
angle of 75°. The asymmetry might be formed by shadowing of the structure by the probe: 
when the probe is placed on the left side of the structure, it reflects the incident light 
and light hardly reaches the full structure; whereas when the probe is positioned on the 
right, the structure is completely illuminated by the beam. Compared to the experimental 
results, the simulations show some differences: Firstly, the wavelengths as well as the z- 
position of the maximum are slightly shifted. This is likely due to the small discrepancy 
between the dielectric constants of the fabricated superlens layers and the data reported 
in the published literature [201121], which are used for the simulations (see Supplementary 
Fig. SI). Secondly, the near-field signal above the topographic step for A > 14.5 fim in the 
simulation shows a saddle shape. The enhanced signal at the topographic steps is possibly 
caused by the sharp edges in the modeled geometry, which give rise to a highly nonlinear 
near- field due to localized corner polariton modes [45j. Please note that this effect is only 
relevant in the propagating polariton regime, whereas it is much less pronounced around 
the superlensing wavelengths since the polariton mode at the flat interface of the step 
becomes also highly localized. In the supplementary information (Supplementary Fig. S3), 
we compare these results with corresponding simulations without the topographic step. We 
find very similar spectral, lateral and vertical response, but the enhanced signal at the edges 
disappears. This comparison verifies that the small topography on the sample interfaces 
due to the sample fabrication has negligible influence on the observed optical signals. The 
enhanced coupling of probe and object by the superlens layer is indeed correlated to the 
material property rather than the topological protrusion. 

We further compare the results of experiment (Fig. 5c) and simulation (Fig. 5d) by 
analyzing the absolute values of the maximum near-field signals with and without SrRuOa 
object, NFw and NF^/^i respectively, as well as the resulting contrast V = {NF^j — 
NF^/^)/{NFw + ^^w/o)- Without object (green curves) we observe maxima at 15.3 /xm 
(c) and 14.5 ^m (d) due to non-localized polariton modes, whereas on the SrRuOa objects 
(red) an enhanced signal is observed for slightly shorter A. The resulting contrasts (yellow) 
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show maxima at 14 fim (c), and 13.3 /xm (d), respectively, which is in good agreement with 
the predicted superlensing wavelength of 13.9 fim. 

Numerical simulations of superlens-enhanced coupling. Why do we observe an 
enhanced signal at a certain distance to the sample surface? The coupled phonon-polariton 
modes, which create the sub-diffraction-limited image, result in a field which decreases 
exponentially with the distance. However, as we place a scattering probe on the image 
side, this probe itself acts as an object next to the superlens creating additional fields at 
the position of the SrRuOs objects. The probe-sample system therefore consists of two 
coupled scatterers with a superlens structure between both of them, showing a resonance 
as a function of the wavelength and the probe-sample distance as discussed in the following 
paragraphs. 

Fig. 6a plots the electric-field distribution of an asymmetric BiFeOa-SrTiOa superlens 
with a line-source in the object plane for different wavelengths and no probe on the image 
side. For wavelengths larger and smaller than the superlensing wavelength around A = 
13.5 /im, we observe unconfined evanescent waves on the image side of the superlens due to 
the excitation of non-localized surface polariton modes at the SrTiOa-air interface. For A = 
13.5 and 14 /im, the superlensing effect takes place with an enhanced confined field on the 
image side of the lens and localized polariton modes at both interfaces of the SrTiOs layer. 
The superlensing effect is further verified by the transfer function simulation of electric 
fields. Figure 6b shows the isothermal contour of the transfer function of the asymmetric 
superlens (for the transfer function of the symmetric superlens see Supplementary Fig. S4). 
Around 13.5 /im, the amplitude of TM waves at the top interface of the superlens still 
maintains reasonably large, even for a tangential wave vector up to lOfco (Fig. 6c). In 
contrast, the field intensity dramatically decays for large wave vectors in a control sample, 
in which the 400 nm SrTiOs film is replaced by a 400 nm BiFeOs layer. 

Figure 6a shows that the electric field at the SrTiOs-air interface possesses in all cases 
a maximum at the sample surface, reflecting the exponentially decaying behavior of the 
polariton mode. The situation changes with the presence of a scattering probe on the 
imaging side of the sample. Fig. 6d shows a simulation for the ideal case of two equally sized 
metal spheres, object and probe, separated by the 2-layer superlens for a fixed wavelength 
of 14 /im and for different gaps z between the top scatterer (i.e., the probe) and the sample 
surface. For small gaps z=25 nm, one observes a large field underneath the probe-scatterer 
Etip due to near- field enhancement as well as a localized polariton mode Eint at the interface 
between SrTiOa and BiFeOa. When the gap is increased to 75 nm, Eup decreases, but Em 
is much larger than that for z=25 nm, whereas both fields decrease with z for gaps larger 
than 100 nm. 

In Fig. 6e, we plot the z-dependence of Eup (green) and Eint (red), as well as of the 
integrated Poynting vector far away from the two objects S (yellow) that corresponds to 
the scattered light intensity measured in our experiments. Eup decreases exponentially 
with the distance, because of the exponential decay of polariton modes at the sample 
surface. However, Em and S show maxima at certain distances between tip and sample 
of 70 nm and 50 nm, respectively. The polariton modes at the interfaces of the superlens 
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enhance the evanescent waves arising from the object. Hence, a larger field amplitude of 
polariton modes Em foi" a certain probe-sample distance indicates enhanced superlensing 
of the coupled probe-object system for this position of the probe scatterer. Moreover, 
the scattered light S from this system has a maximum around the same gap-size, which 
correlates the far-field observation with the localized enhanced coupling. 

Discussion 

The theory of transformation optics [l6l SZl HE] states that a superlens with negative 
index of refraction distorts the optical space in a remarkable way: the space is folded by 
the superlens, with the object plane and the image plane (as well as a plane within the 
lens) at the same position in optical space [46] . If we place scatterers in both planes, one 
can consider them to be two induced dipoles located at exactly the same position in optical 
space, which therefore might lead to the observed enhanced coupling of both scatterers. 

We note that a maximum in the image plane of a superlens was also observed earlier in 
the microwave range using loop antennas as source and detector ^SlISQlEl]. The position 
of the maximum changed due to matched coupling between object and probe, when using 
detector antennas with different radii and different resistance load [50j. However, any 
dependences on the wavelength and on corresponding changes in the sample dielectric 
constants were not studied. It was proposed that this effect could be applicable for three- 
dimensional imaging [50j, which might be possible with the oxide-based superlens as well. 

In this paper, we proposed, designed and demonstrated perovskite-based superlenses for 
electric fields. These materials are in particular suitable for superlenses in the infrared 
range, showing matching pairs of the real parts of their dielectric constants with opposite 
signs. Moreover, perovskites show low intrinsic absorption at the wavelengths of interest 
and different materials can be grown epitaxially on top of each other. Both effects lead 
to low losses in the superlens structure. As some perovskite oxides are ferroelectric (e.g. 
BiFeOa in the current study), a superlens consisting of these materials might be tunable 
by an external electric field. 

We characterized the perovskite superlenses in both spectral and three-dimensional spa- 
tial domain. The near- field examination shows sub- wavelength resolution of A/ 14 at the 
superlensing wavelength. Comparison of symmetric and asymmetric superlens show a 
stronger superlensed signal for the asymmetric case in which the probe is placed directly 
at the superlens interface. This finding is supported by the calculated transfer function of 
the asymmetric superlens, which shows a better performance compared to the symmetric 
case in terms of field enhancement and wave vector bandwidth. 

In addition, we discussed the coupling of near-field probe and object, being placed on 
two opposite sides of a superlens. It is found that the coupling effect is strongly enhanced 
at phonon-polariton resonances and is dependent on the probe-object distance. These 
findings reflect the fact that in a system consisting of object, superlens, and detecting probe 
coupling between all constituents takes place. The observed superlens-mediated interaction 
between two particles might allow for controlled enhanced-coupling effects. These effects 
could find potential applications in local thermal sensors as described by Shen et al [52] 
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as well as in metamaterial-based multifunctional circuits f53l [54] . We envision for example 
an enhanced transport of polaritons from particle to particle by these superlenses: a series 
of superlens-coupled particles can act as a polariton conductor with narrow bandwidth, 
which is determined by the superlensing material. 

Methods 

Growth and absorption losses of perovskite oxides. The various thin films were 
grown via pulsed laser deposition. Thin films of SrRuOa were grown at 680 °C and films 
of BiFeOa and SrTiOa were grown at 700 °C in 100 mTorr of oxygen [55]. Following the 
growth of the SrRuOa layer, features were defined via photolithography and samples were 
ion-milled to produce the objects; subsequently, BiFeOa and SrTiOa films were grown. 
Films were found to be single phase and fully epitaxial in all cases via X-ray diffraction 
and transmission electron microscopy. 

In order to compare the absorption of our perovskite-based superlens with former ma- 
terials used, we calculate the ratios n/ Kj with n and being the real and imaginary parts 
of the refractive index at the wavelengths for which superlensing is expected. Larger n/ n 
corresponds to less material absorption losses. For SrTiOa this ratio is 12.7-10"^, three 
times larger than for SiC and silver, which both show a n/i<i of about 4.2-10"^ at their 
corresponding superlensing wavelengths [18^,^7]. Note that the superlensing condition for 
the perovskite-based superlenses is fulfilled for ssto — —2, whereas for the SiC and silver 
superlens it was observed around e = —3. However, even when esTO = — 3 at A = 14.7 /xm, 
we find nj k — 10.5 • 10~^ being much larger than for the other superlenses. 

s-NSIM setup. In s-NSIM, a scattering probe, namely a metal-coated atomic-force mi- 
croscope (AFM) tip with a typical radius of 50 nm, is placed close to the sample surface 
[34l[35]. This probe transforms the evanescent fields into propagating waves, which can be 
detected in the far-field [34l |56j. In order to separate the near- field from the much larger 
far- field signal, we use the method of higher-harmonic demodulation [56^^ 'STJ: in tapping- 
mode AFM the distance between tip and sample is modulated with small amplitudes of 
about 30 nm at frequencies around Vt — 150 kHz. On this scale, the far-field changes 
linearly with distance resulting in a modulation with the same frequency Vt. On the other 
hand, evanescent waves depend non-linearly on the distance, resulting in a modulation at 
n • (n = 1, 2, 3, ...). Hence, when filtering the modulated signal at n • with n > 2, at 
higher harmonics, we obtain near-field components only, with less far-field contributions 
for larger orders of n. If not stated otherwise, the near-field signals shown in this paper 
represent third-harmonic signals NF^q in order to ensure pure near-field examination. For 
comparisons at different wavelengths, we nomalize NF^q to the current laser power and to 
the spectral response of detector and optical elements in the beam path. 

Our investigation was performed using the free-electron laser FELBE at the Helmholtz- 
Zentrum Dresden- Rossendorf (www.hzdr.de), Germany, which offers continuous tunability 
across a wavelength range of 4 — 250 fim at an average power of up to 10 W (delivered as 
a picosecond pulse train at a repetition rate of 13 MHz). In the 10 to 20 jj^m wavelength 
regime the typical spectral width of the FEL is FWHM=50-100 nm. 
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Near-field interaction and probe-sample coupling. When imaging a superlensed 
signal by s-NSIM, the near-field probe can not be assumed to be a passive element only, but 
the coupling between probe and objects needs to be taken into account. The origin of the 
coupling lies in the optical interaction principle of NSIM, which results in the consequences 
for the NSIM examination of superlenses as discussed in the following. 

In our s-NSIM and the similar one employed in Ref. [18j, the probe-sample system is 
illuminated from the probe side and, hence, the initial near field is generated by the probe 
rather than by the sample. The function of the probe is two-fold: firstly, the evanescent 
waves arising from the probe illuminate the sample locally and excite modes in the sample 
and on its surface. Secondly, the probe senses the fields of the sample modes and transforms 
them by scattering into propagating waves, which can be detected in the far-field. This 
probe-sample interaction depends on the scattering behavior of the probe, the local optical 
properties of the sample, as well as the distance between the probe and the sample surface. 
As a first-order approximation, it can be described by the dipole model introduced by 
Knoll and Keilmann [56j: here, the s-NSIM signal is described by scattering in the near- 
field coupled probe-object system represented through an effective polarizability. A more 
sophisticated approach describes the probe-sample interaction through higher-order modes 
[58j showing the same qualitative results. 

A critical point for successful NSIM measurements is that the probe properties and illu- 
mination are kept constant during the measurements, therefore the changes in the near-field 
signature are essentially caused by local sample properties. In our study we particularly 
keep the position of the probe fixed while the sample is scanned to ensure constant illu- 
mination. The probe properties like material, geometry and orientation with respect to 
the polarization of the incident light are fixed as well. Hence, the observed contrast arises 
from changes in the sample only. The influence of probe or object properties are enhanced, 
when either one of the two is excited close to its resonances. In the mid-IR a metallic probe 
shows no distinct resonances whereas the samples discussed in our work support polariton 
resonances. Hence, it is the sample that determines the NSIM signal, whereas the tip acts 
solely as a small optical dipole. 

In the case of a superlens as the sample the near-field interaction is extended: at the 
superlensing wavelength, the evanescent waves arising from the probe are transferred to the 
object plane. The objects are excited by these waves and generate additional evanescent 
fields, which are reconstructed by the superlens on the image side of the sample. Finally, the 
near-field signal at the image plane is transformed by the probe into detectable propagating 
waves. As the observed signal depends strongly on the properties of the superlens as well 
as the positions of both scatterers, probe and object, with respect to the superlens slab, 
we call this effect a superlens-enhanced probe-object interaction. 

Besides the superlens-enhanced interaction between probe and objects, the probe can 
also interact with the surface layer of the superlens, which is in our case either BiFeOa 
or SrTiOa for the symmetric or asymmetric superlens, respectively. We observe enhanced 
near-field signals due to propagating polariton modes at the toplayer-air interface of the 
samples at wavelengths around 17 /im and 14.6 /im for symmetric and asymmetric su- 
perlens, respectively. Please note that, unlike the localized modes at the superlensing 
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wavelength, this signal does not carry any information about the SrRuOs objects (see 
Figs. 3 and 4). 

Numerical simulations. The numerical simulation results presented in the manuscript 
are all based on the commercial finite-element solver COMSOL 3.5. Due to memory lim- 
itation, the simulation is performed in two dimensions. The extremities of the simula- 
tion domain are assigned scattering properties which essentially mimic the necessary open 
boundary conditions. We modify the size of simulation domain and (local) meshes to ensure 
that the field variation is less than 1 %. Such a convergence of the numerical simulation 
verifies that the boundary condition and meshing are assigned properly. The total mesh 
number is up to 450,000. For a workstation with 32G RAM memory and two dual CPUs 
(2.66 GHz), each simulation takes less than two minutes to converge and reach the relative 
tolerance of 10-6. 

To numerically retrieve the s-NSIM signal shown in Fig. 5, iterative simulations with 
varying tip positions are performed by combining a Matlab script with COMSOL. The 
near-field probe is made by gold, whose dielectric constant is taken by fitting the data from 
Johnson and Christy [59]. The probe is modeled as a triangle, whose tip angle is about 22 
degrees and tip apex is a half sphere with the diameter of 100 nm. The superlens consists of 
layers of BiFeOa and SrTiOa with thicknesses of 200 nm and 400 nm, respectively, and the 
dielectric constants are taken from literature [20l[2T]. The SrRuOa objects are assumed 
to be rectangles of 50 nmx4 fim with a dielectric constant of SrRuOa determined by 
Fourier transform infrared (FTIR) spectroscopy (see Supplementary Fig. SI). The system 
is illuminated by monochromatic light with an incident angle of 75°. In order to study the 
response of the system, we calculate the electric field distribution at all positions of the 
system. We observe a field enhancement underneath the tip and at the superlens interfaces 
due to polariton excitation. The scattering cross sections (Fig. 5b) as measured in the 
experiment are related to the integrated Poynting vector over a solid angle about 70° from 
the tip apex. Such a simulation is repeated for different tip positions, in direct analogy to 
the scanning process. In order to simulate the higher-harmonic demodulation, we calculate 
the gradient of the Poynting vector to obtain the first-harmonic signal. Consequently, the 
gradient of the first harmonic reflects the second-harmonic signal and so on [60j. 
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Figures 




Figure 1. s-NSIM setup and perovskite properties (a) Sketch of the 
experimental setup including the superlens and, the geometry at the near- 
field probe (blue), and the free-electron laser light source. The superlens 
consists of the layers A (BiFeOa) and B (SrTiOa) of thicknesses d and 2d 
((i=200 nm), respectively. The objects to be imaged are SrRuOa patterns 
on a SrTiOa substrate. All constituents of the superlens are perovskite 
oxides that match in their crystalline structure resulting in low scattering 
at the highly crystalline interfaces. The near-field tip probes the evanescent 
fields on the image side of the lens. The superlens is excited by an IR free- 
electron laser, which is precisely tunable in the range from 4 to 250 /xm. 
(b) Imaginary and real parts of the dielectric constants £ of all constituents 
(SrRuOa data determined by FTIR spectroscopy, values for SrTiOa and 
BiFeOa taken from literature jMllSlj, see Supplementary Fig. SI); (c) real 
parts of the dielectric constants at the high-frequency side of their phonon 
resonances depicted in b. The arrow indicates the wavelength at which 
superlensing is expected. 
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Figure 2. Near-field images of three different samples. From left 
to right, the subfigures display: sketches of the sample, topography images 
obtained by atomic force microscopy (scalebars are 10 /im), near- field sig- 
nals as functions of the probe sample distance z for selected wavelengths 
(an offset is added for better comparison), as well as near-field images as 
described in the following, (a) For SrRuOa objects on a SrTiOa substrate 
we image second- and third-harmonic near-field signals {NF2Q and NF^q) 
using a CO2 laser (A = 10.6 /xm). (b), (c) For both types of superlenses we 
depict NFsn at two different wavelengths being A = 17.5 fim and 14.6 fim 
for the symmetric superlens (b) and A = 17.3 ^m and 14.1 /xm for the asym- 
metric lens (c). The red and green curves in the distance curves correspond 
to areas with and without SrRuOa objects on the opposite side of the lens, 
respectively. Such distance curves and the near-field images show strong 
signals at both wavelengths, but only at the shorter wavelengths we observe 
a contrast beyond the diffraction limit due to the superlensing effect. 
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Figure 3. Spectral response of the symmetric superlens. Near-field 
spectrum as well as near-field images and A^Fa^-distance curves for selected 
wavelengths (scalebars are 10 fim). (a) For A = 13.9 to 15.9 /xm, an imaging 
contrast exists due to the localized polariton mode, (b) Shows the near- 
field spectrum for a fixed distance of z=20 nm with the results between 
13.5 /im and 16.25 /im being multiplied by a factor of 4 in the plot, (c) For 
A = 16.8 to 18.4 /im, no imaging contrast is observed although the near- 
field signal is enhanced by the non-localized polariton mode. The red and 
green curves in the distance curves of a,c and the spectrum (b) correspond 
to the near-field signals on areas with and without SrRuOa objects on the 
opposite side of the lens, respectively. 



Figures 



18 




Figure 4. Spectral response of the asymmetric superlens. Near- 
field spectrum as well as near-field images and A^Fa^-distance curves for 
selected wavelengths (scalebars are 10 /im). (a) For A = 13.7 to 15.2 /im, 
an imaging contrast exists due to the localized polariton mode. (b) 
Shows the near-field spectrum for a fixed distance of z=20 nm. (c) For 
A = 15.4 to 18.4 /xm, no imaging contrast is observed although the near- 
field signal is enhanced by the non-localized polariton mode. The red and 
green curves in the distance curves of a,c and the spectrum (b) correspond 
to the near-field signals on areas with and without SrRuOa objects on the 
opposite side of the lens, respectively. 
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Figure 5. Normalized cross sections on an asymmetric superlens. 

(a) depicts experimental data, whereas (b) shows the results from numerical 
simulations (for details see text). The horizontal range in all cross sections is 
12 fim (scalebars are 6 /im). The topography of the sample and the position 
of the SrRuOa object is reflected by the dark areas at the bottom of the 
figures. The localized evanescent flelds on the objects show a maximum at 
a certain distance zq, which increases with smaller wavelength. This effect 
is observed experimentally as well as in the simulations and corresponds to 
a superlens-enhanced coupling of probe and object. (c,d) Spectral behavior 
of the near-field maxima in the experiment (c) and in the simulations (d) for 
areas with (red) and without (green) objects as well as the corresponding 
contrast V (yellow) calculated from this data. The marked areas correspond 
to the observation of non-localized near-field signals (green), and localized 
evanescent fields due to the superlensing effect (red). 
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Figure 6. Simulations of electric field distributions and transfer 
function for the asymmetric superlens. (a) Planar superlens with a 
hne-source as object on one side for different wavelengths, showing a con- 
fined field on the image side of the lens for A = 13.5 — 14 fim due to 
superlensing (all figures with same color scale in arbitrary units), (b) The 
isothermal contour of transfer function in the wavelength range of our inter- 
est plotted versus wavelength A and wavevector kt. The color represents the 
transfer function (the square of the ratio between the transmitted electric 
field after the superlens and the incident field). The white line is the light 
line in air. (c) Transfer functions |rp for the asymmetric superlens (blue) 
and the control sample (red) at 13.5 /im wavelength. The control sample 
replaces the 400 nm SrTiOs film in the superlens by a 400 nm BiFeOa layer. 
One can clearly see that the evanescent wave is enhanced by the superlens 
over a large range of wave vectors (up to lOfeg). The sharp peaks around ko 
are due to total internal reflection, (d) A planar superlens with two spher- 
ical objects on both sides for A = 14 fim and increasing gap z between the 
upper sphere (probe) and the sample surface, (e) Parameters of interest ex- 
tracted from simulations as shown in d: the electric field at the lower apex 
of the probe sphere Eup^ the electric field at the SrTiOa-BiFeOa interface 
Eint, and the integrated Poynting vector S far away from the two-sphere 
system. All parameters are depicted as a function of the gap z. In contrast 
to Etip, which has the highest value at z ^ 0, Eint and S show maxima for 
certain z being 70 nm and 50 nm, respectively. 



